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^ ■ ABSTRACT 

Q"^ ■ We have observed continuum emission at A = 850 /im over ~4 square degrees 

^ ■ of the Ophiuchus star-forming cloud using SCUBA on the JCMT, producing a 

submillimetre continuum map twenty times larger than previous Ophiuchus sur- 
veys. Our sensitivity is 40mJybeam~^, a factor of ~2 less sensitive than earlier 
'nI" . maps. Using an automated identification algorithm, we detect 100 candidate ob- 

\^ ■ jects. Only two new objects are detected outside the boundary of previous maps, 

2 ' despite the much wider area surveyed. We compare the submillimetre continuum 

O ■ Kiap with a map of visual extinction across the Ophiuchus cloud derived using 

^ . a combination of 2MASS and i?-band data. The total mass in submillimetre 

o 

c/3 . mass. A clear association is seen between the locations of detected submillimetre 

cd ■ 

objects and high visual extinction, with no objects detected at Ay < 7 magni- 



> 
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% 



objects is ~ 50 M© compared with pa 2000 M© in observed cloud mass estimated 
from the extinction. The submillimetre objects represent only 2.5% of the cloud 



tudes. Using the extinction map, we estimate pressures within the cloud from 
P/k Ri 2 X 10^cm~'^K in the less-extincted regions to P/k ^ 2 x 10^cm~^K at 
the cloud centre. Given our sensitivities, cold (Td ^ 15 K) clumps supported by 
thermal pressure, had they existed, should have been detected throughout the 
majority of the map. Such objects may not be present at low Ay because they 
may form only where Ay > 15, by some mechanism (e.g., loss of non-thermal 
support). 

Subject headings: stars: formation 



1. Introduction 

Stars must form out of gravitationally bound substructures within a molecular cloud 
but how the substructures themselves form is strongly debated. Possible scenarios for this 
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process range from the gradual release of magnetic support (Mestel & Spitzer 1956; Shu 
1983; Nakano 1984; Shu, Adams, & Lizano 1987; Mouschovias & Ciolek 1999; Basu & Ciolek 
2004), leading to a slow, regulated evolution of structure in the cloud, to violent dissipation 
through interacting waves (Scalo 1985; Mac Low & Klessen 2004 and references therein). 

Rotational transitions of molecules have been important for tracing both substructures 
and motions within molecular clouds. Transitions that probe low densities in clouds (e.g., CO 
(1-0)) reveal widespread emission. Lines observed over low-density regions are themselves 
very wide, indicating the presence of non-thermal motions on large scales that are crucial 
to overall cloud support. Transitions that probe high densities (e.g., NH3 (1,1)), however, 
reveal more localized pockets within clouds, i.e., dense cores. Lines observed in dense cores 
are narrower than those observed in lower-density regions, indicating that non-thermal mo- 
tions are reduced on smaller scales, with a concomitant reduction of support (Barranco & 
Goodman 1998; Goodman et al. 1998). Dense cores are typically associated with embedded 
protostellar objects, suggesting that they are the gravitationally bound substructures from 
which stars form (Tachihara, Mizuno, & Fukui 2000; Di Francesco, Andre, & Myers 2004). 

Although observations of transitions that probe high densities have been a reliable means 
to investigate the small scale substructures of molecular clouds related to star formation, 
observations of submillimetre continuum emission from such regions are an attractive alter- 
native given the high sensitivities such observations have to mass. For example, large-format 
bolometer arrays on large submillimetre or millimetre telescopes have been used in recent 
years to map several pc^ of nearby star-forming molecular clouds, revealing many substruc- 
tures in high-density regions at high sensitivities and resolution (Motte, Andre, & Neri 1998; 
Johnstone & Bally 1999; Johnstone et al. 2000a, 2001). In these regions, the detected con- 
tinuum emission primarily originates from cold dust grains and is almost always optically 
thin, so continuum maps reveal temperature-weighted column densities. These data have 
revealed numerous stellar-mass substructures within the molecular cloud cores whose mass 
distributions have similarities to the stellar initial mass function. 

Most of the regions targeted for widefield submillimetre continuum mapping have been 
those known previously to contain dense cores and active star formation. Much of the 
surrounding regions of lower-density cloud have been ignored in these maps, in part to 
maximize sensitivities. Significant information about the development of substructures may 
exist in these regions, however, but has been overlooked as a result of this bias. For example, 
investigation of Orion B North (Johnstone et al. 2001) found all the substructure to be 
confined to the densest molecular core regions but did not have sufficient cloud coverage to 
provide explicit environmental constraints. To provide an unbiased view of substructures 
within nearby star-forming molecular clouds, we present here submillimetre continuum and 
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infrared extinction maps of ~4 square degrees of the Ophiuchus star forming cloud. These 
observations are part of a larger project, the CO-ordinated Molecular Probe Line, Extinction, 
and Thermal Emission (COMPLETE) project (Goodman 2004), to obtain and compare 
high quality molecular line, submillimetre continuum, and infrared extinction data over the 
extents of the nearby Ophiuchus, Perseus, and Serpens star-forming molecular clouds. (These 
clouds, as well as the Lupus and Chamaeleon clouds and many isolated cores, being observed 
extensively as part of the near- to mid-infrared "From Cores to Disks" (c2d) Legacy survey 
of the Spitzer Space Telescope (Evans et al. 2003).) The data presented here represent 
about half of the entire extent of Ophiuchus that will be mapped in the submillimetre for 
COMPLETE but already they comprise a map that is an order of magnitude larger than 
those of Ophiuchus by Motte, Andre, & Neri (1998) or Johnstone et al. (2000a) (albeit at 
lower sensitivity) and they sample an order of magnitude in column density of the cloud. 

In this Letter, we demonstrate that all substructures detected in the Ophiuchus cloud 
are confined to a small fraction of the submillimetre continuum map. Despite the tenfold 
increase in the area of Ophiuchus probed here, only two objects located off the edges of 
previous maps are detected. No obvious substructures are found below an Ay = 7. Moreover, 
almost all bright submillimetre objects are found in the most-opaque, densest part of the 
cloud where Ay > 15 mag although simple models suggest they would have been detected 
at our sensitivity throughout most of the map if they existed. Instead, a column density 
or extinction threshold for substructure formation, and thus star formation, may exist in 
Ophiuchus and possibly in other molecular clouds. 



2. Observations and Results 

The data were obtained using the Submillimetre Common User Bolometer Array (SCUBA; 
see Holland et al. (1999)) on the James Clerk Maxwell Telescope (JCMT) on Mauna Kea, 
Hawaii^. The area chosen to map consisted of ~4 square degrees in western Ophiuchus where 
Ay > 3 and includes the well-known cores Oph A, Bl, B2, C, and D at its southeastern 
corner. This area was divided up into 38 square fields each ~400 square arcminutes in size. 
Each field was observed over the course of ~1 hour using the "fast scan" observing mode 
of SCUBA where the sky is sampled at a Nyquist rate at 850 /im. Each field was scanned 
six times but at a different chopping amplitude and direction each time. The chop maps 



^The JCMT is operated by the Joint Astronomy Centre in Hilo, Hawaii on behalf of the parent organi- 
zations Particle Physics and Astronomy Research Council in the United Kingdom, the National Research 
Council of Canada and The Netherlands Organization for Scientific Research. 
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were converted into an image by applying the matrix inversion data reduction technique 
(Johnstone et al. 2000b). Note that by using a chopping secondary to map the region, the 
data are insensitive to emission on scales larger than a few times the chop throw, i.e., >120". 
The fields were observed separately through queue mode scheduling throughout Semester 
03A. Since each field was observed at different times and thus under different atmospheric 
opacities, the sensitivity in the final map varies from field to field but only by a factor of 
~1.5. The mean and rms of the 1 a rms sensitivity for all fields is 40 mJy beam^^ and 20 
mJy beam~^ respectively. Since the JCMT beam is ~14" FWHM at 850 /im, the final map 
consists of ~1.9 x 10^ resolution elements, the largest number obtained so far in a single 
SCUBA map. 

Extinctions derived from i?-band star count data by Cambresy (1999) were used for the 
northern half of the map while those derived from Two Micron All Sky Survey (2MASS; 
Skrutskie et al. (1997)) stellar reddening data by Alves et al. (2004) as part of COMPLETE 
were used for the southern half since public release data were available for those fields. Since 
2MASS is an infrared survey, its data are better for probing the highly-extincted regions at 
which the i?-band extinctions saturate. Ay ~ 8. Extinctions derived from i?-band or 2MASS 
data over regions in common show a slight difference of less than 1 mag and a tight linear 
correlation until Ay ~ 8. The resolution of the extinction map depends on the densities of 
stars used to derive extinctions, but it should be only as low as ~3' at the highest extinctions 
(see Lombardi & Alves 2001). The mean and maximum Ay over the surveyed area are 4.2 
and 36.2 mag respectively. 



3. Discussion 

The A = 850 /xm continuum map of Ophiuchus can be used to locate the incidences of 
substructure (i.e., objects of size <120") in the cloud. To locate objects within the map in 
an unbiased manner, the Clumpfind algorithm (Williams, de Geus, & Blitz 1994) was used 
to identify 100 sources of fiux in the map of comparable size to the beam or larger to a 
5 (J threshold in peak brightness. All but 2 of these objects are associated with previously 
identified substructure in the smaller map of Johnstone et al. (2000a) surrounding the known 
Oph cores (i.e., A-G; see also Motte, Andre, & Neri (1998)) Table 1 lists the names, positions, 
and fiuxes of two objects not previously identified since they lie beyond the edges of earlier 
continuum maps. They reside in two new cores in Ophiuchus, which we name Oph H and 
Oph I following Motte, Andre, & Neri (1998). No infrared or X-ray sources toward either 
core are reported within the SIMBAD database. 

Detections of substructure in Ophiuchus allow us to measure their mass and determine 
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their location while the extinction map can be used to estimate the total cloud mass over the 
same area. The integrated flux from all substructure in the cloud is 250 Jy. Assuming a dust 
temperature T^ = 15 K, a dust opacity at 850 /im k^, = 0.02 cm^ g~^, and a distance of 160 pc 
to Ophiuchus then 1 Jy = 0.20 Mq (Johnstone et al. 2000a). Therefore, the amount of mass 
in substructures detected is ~ 50 Mq. From the extinction data and assuming Ayf^n = 2 
X 10^^, the total amount of mass residing within the map area is ~ 2020 Mq. The detected 
substructures, therefore, only account for ~2.5% of the total amount of mass residing in 
the map area. Table 2 shows the distribution of these masses with Ay, and reveals that 
significant mass in substructure is found only at Ay > 15 mag. 

Cumulative mass calculations may be easily biased by a few massive objects, but this is 
not the case for our surveyed area. Figure 1 plots the peak flux, total flux, and radius of each 
identified object versus the Ay at their respective positions. No bright objects are found at 
Ay < 15 despite the large cloud area and mass of that region. A similar result was obtained 
by Johnstone et al. (2001) for submillimetre objects in Orion B North. At high extinction 
the total mass in clumps is determined not by a single object but by an ensemble, although 
the Oph A core accounts for one quarter of the mass in compact objects. Comparison of 
the mass estimation in this paper with the results of Johnstone et al. (2000a) shows that a 
significant amount of mass (almost 50% of the submillimetre fiux in compact substructure) 
in the earlier paper was removed by the fiattening and filtering techniques, stressing the 
difficulty in obtaining proper flux baselines for chopped submillimetre maps. 

What is the physical signiflcance of having most of the detectable substructure only 
at Ay > 15 mag? A threshold in Ay may need to be exceeded for such objects to form. 
Alternatively, our observations may not have been sensitive enough to detect substructure 
in regions of low Ay, e.g. if such objects are larger or less dense than their more-embedded 
brethren. Here, we use a simple model to argue that such objects in locations of low Ay in 
Ophiuchus should have been detected by our observations if they indeed existed. 

Within a molecular cloud, pressure increases with Ay. Following McKee (1989), the 
pressure at depth r is P{r) = 7rG'SS(r), where S is the mean column density through the 
cloud and S(r) is the column density measured inward from the cloud surface to depth r. 
Near the centre of the cloud, S(r) ^ S(s)/2, where S(s) is the column density through the 
cloud at impact parameter s. Thus, pressure in the cloud at depth r can be approximated 
as P{r)/k = 1.7 x 10'^ Ay Ay {s) cm~^ K, where Ay is the mean extinction through the 
cloud, Ay{s) is the extinction through the cloud at impact parameter s, and adopting Ay = 
(S/Sq) mag where Eg = 4.68 x 10^'^ g cm~^. Taking the mean extinction through the cloud as 
Ay ^ 4 and the range of measured Ay = 3 — 30 mag, P/k ~ 0.2 — 2 x 10^ cm~'^ K. Although 
admittedly approximate, the value at the highest extinctions (e.g., the Oph cores) is similar 



-6- 



to the 3 X 10^cm~"^K found by Johnstone et ah (2000a) from fitting stable Bonnor-Ebert 
(BE) sphere models (Bonnor 1956; Ebert 1955) to submillimetre objects in the Ophiuchus 
cores to estimate their masses, internal temperatures, and bounding pressures. 

Pressure support for a molecular cloud must come from non-thermal sources like mag- 
netic fields or turbulence. Most of the previously-detected submillimetre objects in Ophi- 
uchus, however, were fit reasonably well by models of BE spheres supported internally en- 
tirely by thermal pressure (Johnstone et al. 2000a). A similar analysis holds true for Orion B 
North (Johnstone et al. 2001). Note that these fits should not be construed as evidence that 
the objects are in stable equilibrium, since dynamic entities produced in turbulent clouds 
can appear like BE spheres (Ballesteros-Paredes, Klessen, & Vazquez-Semadeni 2003). Sim- 
ilarly, the condensed central regions of clumps evolving via ambipolar diffusion approach a 
similar, albeit flattened, profile (Basu 1997; Basu & Ciolek 2004). The BE sphere model, 
however, provides a lower limit for the critical mass and column density at which an object 
will become unstable to gravitational collapse since only thermal pressure is available to act 
as support. Aside from scaling relations, BE spheres are a one-dimensional family of objects 
defined only by the ratio of central density to surface density, A = p{0)/p{R), with larger 
A implying a greater importance of self-gravity. BE spheres are stable against gravitational 
collapse if A < 14. The mass and density scaling relations of a BE sphere with fixed A are 
Mbe oc P-i/2 and ^BE oc P^/\ 

Can stable objects exist in Ophiuchus where Ay < 15 mag and remain undetected in 
our data? The total submillimetre continuum flux 5*850 associated with a BE sphere is a 
measure of its mass, weighted by temperature. Thus, for BE spheres with a flxed value of A, 
the submillimetre flux as a function of the bounding pressure (assuming that both k,^ and T 
remain constant) is 5*850 oc P(ry^/'^ or 5*850 '^ Avis)^^/"^. Similarly, the peak flux /o measures 
the column density and should scale as /o oc Av{sy^'^ and the radius R oc Avis)^^^"^. The 
total flux, peak flux, and radius increase with increasing A (until a gravitationally unstable 
BE sphere is produced). 

In Figure 1 the above BE sphere relations are overlaid (solid line) on the data to show 
that non-detections at Ay = 10 are incompatible with detections at higher Ay. That is, the 
lack of signiflcant observed substructure at moderate Ay implies that no objects in these 
regions are dominated by self-gravity to the extent of those at higher extinction. Since the 
low Ay regions contain most of the cloud mass, the lack of such detected objects is not due 
to a lack of material. Also shown in Figure 1 is the upper boundary for stable BE spheres 
(dashed-dotted line), assuming that the gas and dust temperature is T = 15 K. Unstable 
objects, those which are capable of collapsing to stars, are only found where Ay > 17. 

How do these results flt with present theoretical models for star formation? If the 
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dominant mechanism for molecular cloud support is magnetic pressure mediated by colli- 
sions between ions (which couple to the magnetic field) and neutrals then the timescale for 
ambipolar diffusion would be long in regions of high ionization, limiting the formation of sig- 
nificant substructure. Crutcher (1999) has shown that the typical magnetic field strength in 
a molecular cloud is strong enough to supply most, if not all, of the required support against 
collapse. McKee (1989) demonstrates that in the lower column density regions of molecular 
clouds, where the extinction to the cloud surface is Ay < 4, the interstellar radiation field is 
capable of maintaining such a high ion fraction. Following this reasoning, the total extinction 
through regions of the cloud where substructure exists should be Ay > 8. This prediction is 
in excellent agreement with our observations. The somewhat higher threshold Ay that we 
find is plausibly a result of the non-uniform structure of the molecular cloud. In the regions 
of high extinction where only cosmic rays are capable of maintaining an ionization fraction, 
ambipolar diffusion enables the growth of prestellar cores. Recent theoretical calculations 
show that MHD turbulence may shorten the timescale for ambipolar diffusion within the 
core by a factor of 3-10 (Fatuzzo & Adams 2002; Zweibel 2002), bringing the ratio of ob- 
served cores with and without stars into agreement with their expected lifetimes. Johnstone 
et al. (2000a) found that approximately half their submillimetre objects in Ophiuchus were 
associated with infrared emission^. 

An alternative model for the support of molecular clouds utilizes turbulent motions 
to supply 'non-thermal' pressure support (Mac Low & Klessen 2004). Such motions are 
observed within molecular clouds via the widths of molecular transitions like CO 1 — 0. In- 
deed, simple virial analysis models indicate that the observed linewidths are large enough to 
provide support of the cloud against collapse, assuming that the 'non-thermal' pressure is 
isotropic and uniform throughout the cloud (Larson 1981). For Ophiuchus typical linewidths 
are AV > 1.5 kms~^ (Loren 1989). Within turbulent clouds, the formation of substructure 
is due to intersecting waves, producing local density maxima, most of which are transient. 
Ballesteros-Paredes, Klessen, & Vazquez- Semadeni (2003) showed that such features can 
look similar to BE spheres despite being non-equilibrium objects. It remains to be deter- 
mined, however, how the formation of substructure depends on local conditions, such as the 
mean pressure. Naively one might expect that these objects should follow a similar scaling 
relation to BE spheres. Regardless of the exact nature of the scaling, it is difficult to under- 
stand how turbulent models can produce the sharp decline in observed objects found in this 
investigation. 

Molecular clouds are likely influenced by both magnetic flelds and turbulent motions. 



^The high sensitivity of the Spitzcr c2d survey wiU hkely change this fraction. 



The SCUBA observations presented in this paper suggest an extinction threshold, which 
is compatible with the predictions of McKee (1989), preventing the formation of relevant 
substructure over the bulk of the cloud. Thus, observations of submillimtere structure in 
molecular clouds yield similar conclusions to surveys of the latter-stage embedded stars (Lada 
& Lada 2003). As with the embedded stars, most if not all of the submillimetre objects reside 
in highly extincted cores, account for about twenty percent of the total core mass, and have 
an IMF-like initial mass function (Motte, Andre, & Neri 1998; Johnstone et al. 2000a). 
Thus, clustered star formation appears limited to the densest regions of molecular clouds 
where turbulent motions and enhanced ambipolar diffusion are able to seed an ensemble of 
prestellar cores. 

The research of D. J. has been supported by an NSERC Discovery Grant. H. K. is supported 
by an NSERC graduate scholarship. This research has made use of the SIMBAD database, 
operated at CDS, Strasbourg, France. 
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Table 1: Newly-Identified Ophiuchus Objects 



R.A. decl. Peak Flux Total Flux Radius Cloud Ay 

Name (J2000) (J2000) (Jy beam-i) (Jy) (") 

H-MMl 16 27 58.3 -24 33 42 OAO 42 36 18 

I-MMl 16 28 57.7 -24 20 48 0.30 2.7 32 7 



Table 2: Percentages of Totals in Ranges of Extinction 
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Fig. 1. — Distribution of the observed properties of measured sub millimetre clumps with 
cloud extinction. Panels from top to bottom show for each clump the integrated flux, peak 
flux, and size. Horizontal dashed lines denote sensitivity and resolution limits while vertical 
dotted lines at Ay = 15 mag mark the boundary of observed bright objects. The solid 
lines show the expected evolution with Ay of a BE sphere with a fixed value of A but only 
marginally observed at Ay = 10 mag (see text). The dash-dotted lines show the upper limit 
for size and mass of stable BE spheres (T = 15 K) as a function of Ay. 



